Storing flux qubits in the quantum RAM of binary Bose-Einstein condensates 
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We consider the hybrid quantum system of an rf-SQUID interacting via its magnetic field with 
a two-component Bose-Einstein condensate. It is shown that the effective coupling between the 
quantized SQUID flux and the atomic hyperfine states scales with the total number of particles 
in the condensate. We demonstrate that this renders the two-component condensate a prospective 
architecture for rapidly accessible storage and retrieval of the quantum information encoded in the 
flux qubits with high transfer fidelity. 
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The concept of hybrid quantum systems, combining 
various advantages (and compensating for deficiencies) of 
either subsystem is rapidly gaining momentum as a con- 
crete tool for the manipulation of quantum information 
Various examples for such hybrid systems include 
cavity-QED with Bose-Einstein condensates (BECs) [||, 
nanomechanical resonators coupled to superconductors 
0, E| j and ultracold atoms coupled to a superconducting 
waveguide cavity 

Using the hybrid system as a quantum computer, one 
part acts as the quantum central processing unit (CPU), 
while the other component is used as the quantum ana- 
log of a random access memory (RAM). A fundamen- 
tal application question is, then, how to quickly store 
quantum information in the RAM which has to possess 
a low decoherence rate, and equally rapidly retrieve it to 
make quantum computations in the CPU. In the present 
study, we investigate the hybrid system shown in Fig.[T] 
and show that it is especially suitable for the particu- 
lar application just described. An rf-SQUID is coupled 
magnetically to a two-component miscible BEC, where 
two hyperfine components, equally populated with N/2 
atoms, interpenetrate. Similar hybrids, integrating cool- 
ing techniques for both ultracold atoms as well as super- 
conductors, have already been demonstrated to be exper- 
imentally feasible 

While the collective Rabi exchange process induced by 
flux qubit operations with the SQUID pseudospin Q flips 
only a single atom pseudospin in the BEC at a time, the 
SQUID-BEC coupling is enhanced by a factor y/N/2 
contributed by each hyperfine component (N is total 
number of atoms in the BEC). This is due to the bosons 
in each hyperfine level being indistinguishable: We do 
not know which of the N/2 single spins has been flipped 
of either component of the BEC. The indistinguishability 
of the atoms in the condensate therefore creates an O(N) 
speedup of the flux qubit transfer process. In addition, 
we will take advantage of the fact that SQUID and BEC 
can rapidly be brought into and out of resonance by tun- 
ing the flux qubit tunneling frequency on very short time 
scales, on the order of nanoseconds [H El- As a re- 
sult, we shall demonstrate that several complete transfer 



cycles of the SQUID qubit state arc possible within the 
decoherence time of the SQUID. The advantages of this 
quantum RAM [H] over other (solid-state based) pro- 
posals for qubit storage in hybrid systems (l34la |. are 
that the quantum information encoded in the hyperfine 
qubit lives for a long time before it decoheres without, 
e.g., the need for further transformation to nuclear spin 
states, and that it can rapidly be transferred back and 
forth between SQUID and BEC. Finally, the quantum 
RAM used in the following has decoherence time scales 
exceeding by several orders of magnitude those currently 
achievable for the parent flux qubits (l9j . 

We first discuss the theoretical description of the two 
subsystems and then derive an effective Hamiltonian, 
taking into account their coupling. The basic assump- 
tion underlying the existence of flux qubits is that the 
flux through a superconducting ring is a quantum vari- 
able governed by an effective potential with two minima 
corresponding to clock- and counterclockwise circulating 
electrons. In the low-energy limit, in which the two-mode 
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FIG. 1. A flux qubit created from superposing left- and right- 
going current states of a rf-SQUID is coupled by its magnetic 
field to the hyperfine transition of a two-component BEC, 
where the Rabi process periodically flips one atom from a 
given hyperfine spin state to the other. The representation is 
not to scale. 
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approximation applies, the SQUID Hamiltonian is |2C 
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where e is the energy difference (the effective flux po- 
tential offset) between the eigenstates |L) and |R) in 
the left- and right-going current basis, and A represents 
the tunneling rate through the potential barrier. In the 
left-right current eigenbasis, a current operator can be 
denned as follows, I ~ Ia z , where I = (L|j($)|L) ~ 
— (R|/($)|R) (exact equality holding when e = 0). The 
current operator is related to the flux operator through 



the SQUID loop by I = 



with $„„ the classical 



externally applied flux and L the loop inductance [211 ]. 
The magnetic vector potential of a loop of radius R in 
the dipole approximation and in polar coordinates reads 
A ~ — ° r a sinfle^,. Thus the qubit magnetic field is 
given by B = V x A — B(r)a z . 

The interacting Hamiltonian for a trapped two- 
component BEC with field operator components ^l/^, 
where the pscudospin a = (\., t) labels the hypcrfine in- 
dex, takes the form (K = 1) 



Hbec = J2 J dDr ^Ur) 
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J * D rH(r)*Ur)*Ar)K(r), ( 2 ) 



where V(r) is the trapping potential, assumed to be inde- 
pendent of spin, Lu a corresponds to the hyperfine energy 
levels, and g a a , are interaction couplings. 

The magnetic coupling of the BEC to the SQUID is 
given by the canonical expression 



(3) 



multiplying the magnetic moment density of atoms with 
the local magnetic field, and integrating over space. Here, 
fi is the total magnetic moment matrix of atoms in the 
a basis. The total Hamiltonian is then 



H = H BEC ® 1 + 1 ® H SQV1D + H b 



(4) 



where 1 is the identity matrix in either BEC or SQUID 
space (corresponding to the first and second vector prod- 
uct entry, respectively). 

We now derive an approximation to the full Hamilto- 
nian, which represents the coupling of the flux qubit to 
a large two-component BEC on the mean-field level. To 
this end, we first truncate the field operators using the 
standard mean-field approximation to just one mode for 
each species a, 



(5) 



where 4> a (r) are mean- field solutions for the interacting 
BEC, leading to coupled Gross-Pitaevskh equations. 

As we will see below, even with bosonic enhancement, 
the SQUID-BEC coupling is small compared to the hy- 
pcrfine energy difference, so that their ratio is a small 
parameter. In terms of this parameter, the two systems 
are still weakly coupled. Along with the form of the inter- 
action ([3]) , which only couples states in the BEC Hilbert 
space that differ by a single pseudospin, this leads us to 
conclude that the full Hilbert space of the BEC can be ap- 
proximated by a two-dimensional subspace; the "ground 
state" with and atoms and an "excited state" that 
differs by only a single atom of each spin. Wc therefore 
restrict ourselves to the following two BEC basis states, 



(4P(4)> ac > 




V(AW)V(^t + l)! V 1 



(6) 



In the weak coupling limit (in the sense described above), 
the dominant process is the periodic transfer of a sin- 
gle quantum of energy from one subsystem to the other, 
i.e., the Rabi cycle. The whole binary BEC appears in 
this limit, and within mean- field Eq. (|5|), as a giant two- 
level atom coupled to one "cavity" mode produced by 
the SQUID. We have, in addition, verified explicitly that 
processes flipping more than one spin at a time have neg- 
ligibly small probability. 

As the BEC Hamiltonian ([2]) conserves spin, its 
off-diagonal elements in the basis of ((6]) vanish, 
while the diagonal ones give the mean-field energies 
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Using that £ Nf+1 



A' i 



lu^ — lo^. (the interaction renormalization of the hypcr- 
fine splitting being small compared to the level spacing), 
the interacting Hamiltonian ((4]), can finally be written 
as, omitting 0(1/ N) and constant terms, 
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The SQUID-BEC coupling matrix is expressed as 
9a,,' = J d D rfi(r)B(r)i(r), 



(8) 



and e' — e + N^g^. ^ ■ A*t,t + ^iffj. i ' 1*1,1 i s a renormalized 
level splitting of the SQUID. 

A first major result of the present investigation is the 
0(N) enhancement of the effective coupling in the last 
term of Eq. 0. Using this effective Hamiltonian, it is 
easily shown that one can prepare the flux qubit in an 
arbitrary state and then pass it to the BEC; because of 
Rabi oscillations this state will then get passed back to 
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the SQUID, with the process repcatable for many cycles 
due to the short duration of transfer. We first focus on 
controllcdly passing the state to the BEC for it to remain 
there. This can be achieved by bringing the SQUID into 
and out of resonance with the given hyperfme splitting 
of the atoms. Here, we will crucially take advantage of 
the fact that one can dynamically change the flux tun- 
neling barrier of the SQUID potential on a essentially 
arbitrarily rap id timescale within a large frequency win- 
dow (Ref. [10( quotes a change of several GHz for A on 
timcscales of nanoseconds). Changing the tunneling rate 
amounts to replacing A — > A(i). The time dependence 
should be chosen such that the SQUID and the BEC are 



only in resonance, i.e., A(i) 



lo^ for exactly half 



a Rabi period when one transfers the flux qubit to the 
BEC. Outside of this window they should be far off res- 
onance. One choice for this is simply 



A(t) = A W(t), 



(9) 



where A n 



lj^ — lo^ and W(t) ramps from an arbitrarily 



chosen off resonance value to 1 for half a Rabi period and 
then ramps back to its off resonance value. To remain in 
the qubit subspace, i.e., to not excite higher SQUID flux 
potential levels, the ramping time should remain slow 
compared to A^ 1 . 

Assuming the SQUID is at the degeneracy point for 
simplicity, e' = 0, the energy cigenstates of an isolated 
SQUID are |0) s = ^[|L)+|R>] and |1} S = ^[|L)-|R>]. 
The computation basis for the combined system is then 
taken to be |00) = |0) B <8 |0) S) |01> = |0} B 8> |l)s, 1 10) = 
|1)b®|0)s, and |11) = |1)b®|1)s- For a system with equal 
intial populations, 2V1 = N, = N/2, the Hamiltonian in 
the computation basis reads, setting A = w. 
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where Q = ^9 ' Mt4 * s enhanced by a factor N/2 over 
the single-particle coupling. 

The Hamiltonian (fTUfl is easily transformed into a 
cavity-QED type Hamiltonian. Here the two-component 
BEC corresponds to a two-level atom, while the SQUID is 
the cavity photon. For the SQUID we let — a z — > Aa^a, 
where a, a' are photon annihilation and creation oper- 
ators. Assuming a real Rabi frequency Q, we replace 
<r x — > a' + a for the SQUID and a x — > a + + a_ for the 
BEC operators in the coupling term. We then have in 
rotating wave approximation [22| 



We take as a concrete example a 87 Rb BEC 
of N — 10 6 atoms, which has been shown to 
give miscible two-component condensates when the 
states \\) = \h 2 Si/ 2l F = l,m F — -l) and |t) = 
\5 2 S 1/2 ,F = 2,m F = -2) are used 0. For the pur- 
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pose of our computation examples below, we set A 
6.835 GHz equal to the hypcrfinc splitting without Zee- 
man shift between \l) and |t) (a non-degenerate two-level 
hyperfme qubit can be formed from these states by ap- 
plying a small external magnetic field to lift the rap de- 
generacy) . The single-particle coupling strength is taken 
to be \g ■ ~ 20 Hz, corresponding to a SQUID ra- 
dius of 1 fxm, with a current of 1 fiA, and a BEC-SQUID 
distance of about 20 /im. The fact that the latter dis- 
tance is much larger than the SQUID radius implies that 
the detrimental effect of Meissner screening on the trap- 
ping is negligible 0] , representing a major advantage over 
superconducting cavities, where the cold atom cloud is 
situated close to a large waveguide geometry The 
spin-flip lifetime due to current fluctuations should not 
be a limiting factor at distances of this order as well [25| ; 
finally, the SQUID needs to be cooled to temperatures 
below 100 inK, to avoid disturbance of the qubit genera- 
tion and transfer process by thermal photon emission. 

An arbitrary state in the computation basis can be 
decomposed as the superposition 



I*) 



Cl |ll) 



c 2 |01) + c 3 |10) + c 4 |00}. (12) 



The coefficients c i satisfy the coupled Schrodingcr equa- 
tions d t Ci{t) = —iH(t)ci(t), for a given Cj(t = 0) initial 
condition. We now provide two examples for a trans- 
fer protocol, where the results are calculated numerically 
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H = -^■cr 2 ®l + l(X)Aa t a-fi(cr_®a t +(T + ®a), (11) 

which is a Jaynes-Cummings Hamiltonian cf., e.g., [23| . 
For the calculations presented below, we have verified 
that (a' a) < 1 at any instant, so that at most one photon 
is in the SQUID mode. 



FIG. 2. Top panel: Probabilities for finding the system in 
the state |01) and 1 10} as a function of time for two different 
ramping times. Left: Ramping time 5 ns with final (time- 
averaged) fidelity F = 99.3%. Right: 10 ns ramping with 
final F = 95.9% (we fix A^ 1 = 0.15ns). Bottom panel: Final 
time-averaged fidelity as a function of ramp time. 
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FIG. 3. The top panel shows the transfer of a superposition 
state of the flux qubit \ip(t = 0)) = |0) B ® 2~ 1/2 (|0) s + |1> S ) = 
2~ 1/2 (|00) + |01>) to the hyperfine qubit state \ip(t f )) = 
2" 1/2 (|0) B + |1>b) <8 |0)g = 2~ 1/2 (|00) + |10». The middle 
panel shows the production of the state \tji) = 2 _1//2 (|01) + 
1 10)) from an initial state |01). Finally, the bottom panel 
demonstrates a memory transfer process with a final fidelity 
of 99%. The ramping times for all figures are 5 ns. 



without using the rotating wave approximation. 

We start with a simple transfer protocol. If at time 
t = the combined system is in the state \^>) = |01), then 
c 4 (f = 0) = for i ^ 2 and c^f = 0) = 1. The probability 
that at some later time the system is found in the state 
1 10} is c 3 (i)| 2 , see Fig.[2j where we give examples of the 
fidelity F = K^f/ (target) |<F)|, evaluating for this protocol 
to F(t) = \aQ (t)c 00 (t) + aio(t)c 10 (t) | writing the time 
evolution of the initial state as |\E') = ^ m n c m,n{t)\ mn ) 
and similarly that of the target state as |>F (target)) = 
S m n a m n WI TOn )- We define ramp time as the time in- 
terval it takes to raise A from 1 % above its value at 
t — > — oo to within 1 % of being on resonance, using a 
smooth tanh profile for W(t) [26j. The bottom panel of 
Fig. [5] shows that F decreases with the ramp time thus 
defined. Rapid transfer (while staying within the limits 
of the two-mode approximation) therefore significantly 
improves on the fidelity. 

More elaborate examples for flux qubit storage and re- 
trieval are displayed in Fig. [31 where we show in the top 
panel that an almost perfect transfer (fidelity F = 99.3%) 
of a flux qubit to the hyperfine qubit is possible within 
a short time (we have verified that an initial superpo- 
sition state with phase difference 7r/2, \ip(t = 0)) = 
|0) B ®2- 1 / 2 (|0) g +i|l) s ) = 2- 1 / 2 (|00) + i|01)) yields the 
same result). The middle panel of Fig. [3] shows that one 



can entangle the two qubits, i.e., produce a final state 
= 2~ 1 / 2 (|01) + 1 10) ) , by preparing the system in the 
state |01) at t = 0, and then bring them into resonance for 
a quarter of a Rabi cycle. The latter process is in the ex- 
ample shown achieved with a fidelity of 98.8 %. Finally, 
the bottom panel demonstrates the quantum RAM in 
storage and retrieval action, by transferring first the ini- 
tial qubit state to the BEC and then back to the SQUID 
by again tuning into resonance, showing that this con- 
crete RAM-type process can be achieved as well with 
large fidelity, here F = 99%. 

In conclusion, we have demonstrated that rapid and ef- 
ficient flux qubit storage and retrieval are feasible when 
an rf-SQUID is magnetically coupled to the hyperfine 
levels of a two-component BEC. We expect a periodic 
extension of the present hybrid to be conceivable, with 
atomic ensembles in a large-filling optical lattice located 
on top of an array of SQUIDs. With such a setup, the 
flux qubits can be entangled with each other via the con- 
densates in the lattice sites. 

Transferring back and forth flux qubits to long-lived 
hyperfine qubit states of atoms in a two-component BEC, 
for which one has rather precise control over the micro- 
scopic quantum physics, illuminates further whether the 
postulated quantum variable flux has fundamental char- 
acter. This is directly related to long-standing intriguing 
questions on the existence of macroscopic quantum vari- 
ables and thus the quantum-classical transition 27, 28| . 
Among future research directions, then, is the potential 
influence of microscopic variables of the SQUID on the 
qubit transfer process. 
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